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Abstract. The coherent physical alignment of galaxies is an important systematic for grav- 
itational lensing studies as well as a probe of the physical mechanisms involved in galaxy 
formation and evolution. We develop a formalism for treating this intrinsic alignment (lA) 
in the context of galaxy-galaxy lensing and present an improved method for measuring lA 
contamination, which can arise when sources physically associated with the lens are placed 
behind the lens due to photometric redshift scatter. We apply the technique to recent Sloan 
Digital Sky Survey (SDSS) measurements of Luminous Red Galaxy lenses and typical (~ L*) 
source galaxies with photometric redshifts selected from the SDSS imaging data. Compared 
to previous measurements, this method has the advantage of being fully self-consistent in 
its treatment of the I A and lensing signals, solving for the two simultaneously. We find an 
lA signal consistent with zero, placing tight constraints on both the magnitude of the lA 
effect and its potential contamination to the lensing signal. While these constraints depend 
on source selection and redshift quality, the method can be applied to any measurement 
that uses photometric redshifts. We obtain a model-independent upper-limit of roughly 10% 
lA contamination for projected separations of ~ 0.1-10 /i~^Mpc. With more stringent 
photo-z cuts and reasonable assumptions about the physics of intrinsic alignments, this up- 
per limit is reduced to 1-2%. These limits are well below the statistical error of the current 
lensing measurements. Our results suggest that lA will not present intractable challenges to 
the next generation of galaxy-galaxy lensing experiments, and the methods presented here 
should continue to aid in our understanding of alignment processes and in the removal of lA 
from the lensing signal. 
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1 Introduction 

Observations of gravitational lensing, the deflection of light by matter between the source and 
observer, have become an important and widely-used tool in cosmology and astrophysics (for 
a review, see [1-4]). Lensing measurements are equally sensitive to all types of matter, making 
them a particularly useful probe of dark matter and theories of gravitation [5]. Galaxy-galaxy 
lensing measures the gravitational influence of massive foreground ("lens") objects on the 
measured shapes of background ("source") galaxies. In the weak lensing regime, such effects 
must be studied statistically, since the lensing deflections are typically factors of several tens 
to hundreds smaller than the intrinsic galaxy ellipticities. Correlating the observed source 
shapes with lens positions measures the galaxy-mass cross-correlation, which can in turn 
be used to examine the density profile of halos and probe the standard ACDM paradigm 
as well as theories of modified gravity [6, 7]. Combining such lensing results with clustering 
measurements of similar lenses may provide an especially robust probe of density fluctuations 
[8]. 
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Coherent, large-scale correlations of the intrinsic shapes and orientations of observed 
galaxies are a potentially significant source of systematic error in gravitational lensing stud- 
ies, with the corresponding potential to bias or degrade lensing science results. This intrinsic 
alignment (lA) has been examined from both observational [9-14] and theoretical [15-27] 
perspectives. For instance, luminous red galaxies have been observed to align preferentially 
towards overdense regions. Since lensing increases the observed tangential shape of back- 
ground galaxies, contamination from such radial alignment would reduce the inferred lensing 
effect. This type of lA contribution is often referred to as a "GI" term, referring to the 
correlation between matter density, which sources gravitational (G) shear, and intrinsic (I) 
shear (see [23, 26]). It is this type of lA effect that can contaminate a galaxy-galaxy lensing 
measurement, and one can relate the observed distribution of lens galaxies to the underly- 
ing matter distribution to equate with the GI term as originally formulated in the context 
of cosmic shear. Thus, with sufficient modeling of lA and the effects of redshift quality, a 
measurement of lA in galaxy-galaxy lensing can probe the lA contribution to cosmic shear. 

In the context of galaxy-galaxy lensing, the use of accurate spectroscopic redshifts would 
ensure that lens and source galaxies are separated by a large distance along the line-of- 
sight and thus that any intrinsic correlations are negligible. Although such correlations 
may exist between the shapes of the background source galaxies, they yield no total signal 
when measurements around many lenses are "stacked." However, obtaining large numbers of 
spectroscopic redshifts is highly resource intensive, making it difficult to form large samples 
and achieve small statistical uncertainty. An alternative approach is to use photometric 
redshifts (photo-z's) from multi-band imaging. Although less accurate than spectroscopic 
redshifts, photo-z's make much larger galaxy catalogs readily available, greatly improving the 
potential signal-to-noise of lensing measurements. The degradation of distance information 
due to the use of photo-z's requires additional calibration during lensing analysis [28, 29] 
and may introduce additional systematic uncertainties. As many of the next generation 
surveys will be photometric (e.g., KIDS^, DES^, HSC"^, and LSST^), a detailed quantitative 
understanding of the impact of using photo-z's will be crucial for the future of precision weak 
lensing. 

When using photo-z's in a galaxy-galaxy lensing measurement, lA can contribute a 
non-negligible signal. Large redshift uncertainties can result in galaxies that are physically 
associated with the lens being mistakenly assigned a large line-of-sight separation. Thus, if 
there exists a significant correlation between galaxy shape and the density field, as traced by 
the location of lens galaxies, lA contamination can enter the lensing signal. 

Despite previous work, the physical mechanisms underlying lA, as well as the magnitude 
of lA in relevant lensing samples, remain poorly understood. In this paper, we develop and 
apply a technique to measure this lA contamination. Previous studies have used low-redshift 
spectroscopic observations to directly measure the lA shear - density correlation (e.g. [10, 30, 
31]). In [13], a similar measurement is made using photometric redshifts of Luminous Red 
Galaxies (LRGs) at intermediate redshift (0.4 < z < 0.7), requiring the removal of the lensing 
effects that contribute due to redshift uncertainty. In [32], lA is constrained in the context 
of galaxy-galaxy lensing under the assumption that no lA contamination exists in a photo-z 
selected group of background source galaxies. We propose a new method to simultaneously 

^KIlo Degree Survey, http://kids.strw.leidenuniv.nl/ 
■^Dark Energy Survey, littps://www. darkenergysurvey.org 
^Hyper Suprime-Cam, http://oir.asiaa.sinica.edu.tw/hsc.plip 
''Large Synoptic Survey Telescope, http://www.lsst.org 
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measure the lA and lensing signals from photometric lensing measurements with minimal 
assumptions. We exploit the fact that the lA signal is sourced by contamination from galaxies 
physically associated with the lens. Splitting the source catalog by separation in photo-z from 
the lens allows us to compare samples with different levels of contamination and thus solve 
simultaneously for the lA signal and lensing shear. This method is completely self-contained 
and model-independent, using only the lensing measurements themselves to constrain lA. 
Limited model-dependent assumptions can also be applied to improve constraints (see [33- 
36] for techniques to remove lA using flexible models). 

We apply this method to LRG lenses and typical photometric galaxy lenses in the Sloan 
Digital Sky Survey (SDSS). The main distinctions between this study and [13] are the lower 
redshift range for this work, and the fact that [13] directly measured the intrinsic alignments 
of LRGs whereas we use the galaxy-galaxy lensing signal to measure the intrinsic alignments 
of typical source galaxies that dominate lensing measurements. Using this technique with 
upcoming large imaging surveys will allow us to probe lA across a range of scales and as a 
function of both lens and source properties. A potentially valuable application of this method 
is in the galaxy-galaxy lensing mass determination of galaxy clusters. The use of clusters to 
constrain cosmological parameters relies on understanding and removal of effects which could 
bias the inferred masses. Better characterization of alignment in different galaxy populations 
will improve our understanding of systematic errors for weak lensing studies in general, and 
will provide a probe of the complex astrophysical processes involved in galaxy formation and 
evolution. 

The paper is organized as follows. Section 2 describes the observations involved in 
galaxy-galaxy lensing and the data utilized in this study. Section 3 contains a brief overview 
of the formalism of galaxy-galaxy weak lensing, including lA. Section 4 discusses our method 
of extracting lA from the lensing signal, and section 5 shows our results. We conclude with a 
discussion of the implications of our work in section 6. An appendix includes a more technical 
discussion of issues relevant to our measurement method. Calculating the lensing signal and 
calibration factors requires a cosmological model, although this assumption has a minimal 
impact on the lA results. Where relevant, we adopt a flat ACDM universe with Qrn = 0.25 
and present results in units with h = 1. 

2 Data 

For measurement of galaxy-galaxy lensing and of intrinsic alignments, we use data from 
the Sloan Digital Sky Survey. A requirement for this study is to have reasonably accurate 
redshift information for the galaxies that we use as lenses (tracing the density field) and 
those that we use to measure shapes (tracing either the intrinsic alignments or the lensing 
shears, depending on the line-of-sight position with respect to the lens). In addition, we must 
measure galaxy shapes at sufficient resolution. In this section, we describe the data used for 
the lens and shape samples. 

The SDSS [37] imaged roughly vr steradians of the sky and spectroscopically followed-up 
approximately one million of the detected galaxies [38-40]. The imaging was carried out by 
drift-scanning the sky in photometric conditions [41, 42], in five bands (ugriz) [43, 44] using a 
specially-designed wide- field camera [45] . These imaging data were used to create the source 
catalog used in this paper. We also use SDSS spectroscopy for the lens galaxies. All of the 
data were processed by completely automated pipelines that detect and measure photometric 
properties of objects, and astrometrically calibrate the data [46-48]. The SDSS I/II imaging 
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surveys were completed with a seventh data release [DR7, 49], though this work also relies 
on an improved data reduction pipeline that was part of the eighth data release, from SDSS 
III [50]; and an improved photometric calibration ['ubercalibration', 51]. 

2.1 Lens sample 

The lens sample considered here consists of LRG lenses from DR7 as selected by [52]. To 
match cuts in the shape sample on imaging and PSF quality, positions with respect to bright 
stars, and extinction, we remove 8% of the area, yielding a final area of 7,131 deg^. The 
lens sample consists of 62,081 galaxies with spectroscopic redshifts between 0.16 < z < 0.36 
and (7-band absolute magnitudes Mg in the range [—23.2,-21.2], with a roughly constant 
comoving density of 10~^(/i~^Mpc)~'^. For all calculations, we include a weight for each lens 
designed to mitigate the effects of fiber collisions, completeness, and large-scale structure 
fluctuations [52]. 

2.2 Shape sample 

The sample of galaxies with shape measurements, used in this paper to study both the 
gravitational shear and the intrinsic alignment field, is described in [53]. This catalog is based 
on the SDSS DR8 photometric data, processed using re-Gaussianization [54] to correct for 
the effect of the point-spread function (PSF). It contains ~ 30 million galaxies (1.2 arcmin"^) 
with r-band apparent magnitude r < 21.8. There are also cuts on the image quality, quality 
of the PSF estimation, galaxy size compared to the PSF, and Galactic extinction. The 
photometric redshifts assigned to each galaxy based on the five-band photometry are from 
the Zurich Extragalactic Bayesian Redshift Analyzer [ZEBRA, 55]. The photo-z uncertainty 
is iTz/ (1+z) ~ 0.11, due primarily to the low S/N for the majority of the photometric galaxies. 
The impact of this uncertainty on lensing measurements with this catalog is characterized 
by [29]. For this analysis, we use the entire source sample, as well as "red" and "blue" color 
sub-samples, where the split is done using the ZEBRA template type (see appendix A for 
more details). For comparison with other I A results, note that the typical luminosity of these 
sources in the redshift range of interest is ~ L*, making them comparable to the SDSS L4 
sample. Where relevant, we refer to the SDSS luminosity classification scheme described in 
[56], based on r-band absolute magnitude. The labels {LI, L2, L3, L4} correspond to Mr in 
the ranges {[-17,-16]; [-18,-17]; [-19,-18]; [-20,-19]}. 

3 Lensing formalism 

As light travels from a distant galaxy, the presence of intervening matter alters its path, acting 
as a "gravitational lens" and altering the observed shape of the galaxy. This gravitational 
shear is determined by the projected mass density along the line-of-sight. In this section we 
summarize the relevant aspects of lensing formalism and develop a consistent treatment of 
intrinsic alignment. For a more detailed treatment of lensing, see, e.g., [1] or [53, 57] for 
systematic issues. 

3.1 Galaxy- galaxy lensing 

Observed galaxy shapes are a combination of the intrinsic shape and the gravitational shear. 
The observed ellipticity of an object, e, is related to observed shear by the responsivity 
factor TZ, which measures the average response of an ensemble of source shapes to a given 
shear: (7) = {e) /2TZ [58]. In the weak lensing limit, one can write the shear as a sum of 
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gravitational (G) and intrinsic (I) contributions: 7 = 7*^ + 7^. For an individual object, 
the intrinsic shape dominates, although if the intrinsic components of galaxy shapes are not 
correlated with each other, they will average to zero and simply contribute a "shape noise" to 
the measurement. Thus, it is often assumed that the observed shear components provide an 
unbiased estimator of the true gravitational shear: over a large number of sources (7) = 7*^, 
where a tilde denotes an observed quantity. In the case of galaxy-galaxy lensing, the lensing of 
background ("source") galaxies is dominated by a single massive foreground ("lens") object, 
which imparts a tangential shear to the observed shapes, yielding a correlation with the lens 
position. We define the lens surface density contrast, AS, as the difference between the 
average (projected) surface mass density within radius Vp and the surface mass density at r^: 
AS(rp) = S(< Vp) — S(rp). The tangential gravitational shear, 7^^, is related to the surface 
density contrast: 

AS = Se7f , (3.1) 
where Sc is the critical surface density, expressed in comoving coordinates as 

cP D, 



^T^G {I + zif DiDi, 



(3.2) 



The angular diameter distances Di, Ds, and Dig are between the observer and lens, observer 
and source, and lens and source, respectively, while zi is the lens redshift. 

If the intrinsic shapes of nearby galaxies are correlated with each other, an intrinsic 
alignment shear contribution, which we denote 7^^, will be present: 

(7t) = 7? + 7^^ (3.3) 

For notational simplicity, in the remainder of this work we drop the "t" subscript - all shears 
are assumed to refer to the tangential component. A galaxy- galaxy lensing measurement 
consists of averaging over many lens-source pairs to estimate the average lens density contrast: 



lens lens 



lA 



B{rp) E Wjtcj^j B{rp) Wjtcj (^AS(rp)Scj ^ + 7^ 

^(^^) = = ' (3.4) 

E E Wj 

j j 

where wj is the weight given to each lens-source pair, and B[rp) is the boost factor, dis- 
cussed below, which accounts for sources that are physically associated with the lens. In the 
summation expressions, "lens" denotes a sum over all real lens-source pairs with projected 
separation in the desired bin, while "rand" denotes a sum over random lens-source pairs 
(where the random lenses are distributed with the same angular and radial selection function 
as the real lenses). The optimal weight [28] for each pair is a combination of the geometric 
factor Scj-, shape noise from the variance in the source ensemble (erms)) and the individual 
object measurement noise {(Tej)- 

Wj = y — ^. (3.5) 

c,j y^ims I ^e,j 
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3.2 Accounting for physically associated galaxies 

When looking along the line-of-sight at or near a lens galaxy, the measured number density 
of sources, as a function of redshift, is given by the sum of two components: a smooth 
background dn/dz, determined by observational parameters and selection cuts; and a sharp 
peak located at zi due to galaxies that are physically associated with the lens and thus 
strongly clustered with it. The boost factor, B, measures the contribution of this excess peak 
in the photometrically defined source sample. Since these excess galaxies are not lensed, they 
will dilute the measurement unless accounted for explicitly, as done in eq. 3.4. The boost 
is observationally determined by comparing the weighted number of lens-source pairs for a 
given source sample with the number of random-source pairs for the same sample: 

lens 

^(^^) = 7h— (3.6) 

j 

With accurate redshift information, it would be easy to exclude physically associated sources, 
and the boost would approach 1. The measured boosts for the samples used in this study 
are shown in figure 1. The quantity {B — 1)/B gives the fraction of galaxies in the source 
sample that are physically associated with the lens. Note that we are assuming that the 
number of physically associated galaxies is the same as the number of "excess" galaxies due 
to clustering. This assumption need not hold, a subtlety we will explore below. 

3.3 Correcting for photometric redshift bias 

If lA contributes no signal, = 0, and the expression in eq. 3.4 can be rearranged: 

lens _ rand 

^^IA=0 _ j _ j = A I 1 f'i 7\ 

AS ~ lens ~ rand - + ^, {'i-l) 

E E U)j 

3 3 

where we have defined b^, the bias in the lensing signal due to the uncertainty in photo- 
metric redshift, and used the fact that excess galaxies are at the lens redshift and thus have 
Scj~^ ~ such that 



lens rand 

E 

j 



Because the lensing efficiency, expressed by Sc, is a nonlinear function of distance, even an 
unbiased scatter in redshift can result in a biased lensing measurement, and correction for 
this photometric bias must be applied when photometric redshifts are used. The value of 
bz is a function of the lens redshift, but for the purposes of this study, we only need the 
average value (bz), determined by integrating over the lens redshift distribution (see, e.g., 
eq. 22 of [29]). Combined with the definition of the boost, we can express bz as a sum over 
sources around random points rather than lenses, allowing the use of calibration sets from 
fixed areas of the sky unassociated with the lenses. To calculate (bz), we follow the methods 
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Figure 1. Boost factors, which measure the number of "excess" lens source pairs due to clustering 
and photo-z uncertainty, are shown for different source colors and source subsamples (described in 
sec. 4). For red sources, the "extended boosts" (described in sec. 4.2) are shown as open symbols. 
As expected, boosts are larger for red galaxies (which cluster more strongly) and for samples defined 
with photo-z's closer to the lens redshift. In all figures, error bars indicate 68% confidence regions 
and small horizontal offsets are added for clarity. 



of [29], utilizing the same calibration sets. These calibration sets consist of galaxies with both 
spectroscopic and photometric redshifts, allowing the determination of both T,c,j and '^cj- 
Object selection in these calibration samples matches that for the source samples considered 
here. 

With these definitions, we can rearrange eq. 3.4 to write the true AS in terms of the 
observed AS and the lA contamination (no longer assumed to be zero): 



AS( 



AS(r„ 



lens \ 

i 

lens 

j 



(3.9) 



where = ((5^) + 1) . In the limiting case of no lA, physically associated galaxies con- 
tribute nothing to the measured signal, and the boost acts to cancel this "dilution." 
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4 Methodology for measuring intrinsic alignment 



4.1 Solving for I A 

Galaxies that are physically associated with the lens but mistakenly placed behind the lens 
due to scatter in photometric redshift can lead to lA contamination. At greater line-of-sight 
separation from the lens, the number of associated galaxies should decrease, assuming a well- 
behaved photo-2: distribution. We exploit this fact to separate lA from the lensing signal. 
Our I A method is valid in the presence of catastrophic photo-z failures, however the rate and 
distribution of these failures must be understood using a representative redshift calibration 
sample. In principle, catastrophic failures could induce a non-physical correlation between 
I A and observed line-of-sight separation which will bias results if not taken into account. For 
these SDSS measurements, the failure rate is both low and well-understood and can be safely 
ignored. 

With a minor simplification to the second term in brackets in eq. 3.9, we can explic- 
itly solve for the lA contamination. As discussed above, lens-source pairs can be divided 
statistically into "random" and "excess" pairs: 



lens 



E 



rand 



(4.1) 



We take the sum over "random" pairs to be zero, since non-associated pairs should not be 
aligned with the lens, and it is assumed that the number of excess pairs is a good approx- 
imation for the number of physically associated ones. On small scales, where clustering is 
strong, this assumption is a good one. On larger scales, this assumption may break down, a 
possibility we address below. We further simplify the expression by replacing the value 
for each lens-source pair in the sum with the effective average value for excess pair, 7ia('"p): 



lens 



(4.2) 



Equation 3.9 can then be rewritten: 



AS(rp) 



AS(rp)-(i?(rp)-l)7iA(rp)(Se)c 



(4.3) 



where {T,c)ex{rp) can be directly measured: 



(Ec)ex(^p) 



ixcess _ all _ rand 

_2 J 3_ 

excess ^ rand 



(4.4) 



This simplification of the lA sum is reasonable because the large photometric uncertainty 
effectively removes any correlation between observed Sc and for each physically associated 
object. In other words, the intrinsic alignment shear should depend on the line-of-sight 
separation between lens and source, but our photo-z errors are far larger than the scales 
relevant for lA variation. 
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We consider two groups of lens-source pairs, denoted a and b, defined in terms of the 
separation between lens (located at zi) and source (with photometric redshift Zp). We take 
sample a to consist of all pairs with zi < Zp < zi + Az while sample b has Zp > zi + Az, 
where Az = 0.17. Other splitting schemes are feasible, and this technique can be generalized 
to involve more than two source subsamples, although given the limitations in the statistical 
power of current measurements, using two subsamples is preferred. The potential advantage 
of using additional subsamples is discussed in section 6. We assume here that the lens redshift, 
zi, is determined precisely with spectroscopic measurements, but it is straightforward to 
include uncertainty in zi. We also refer to the src and assoc samples, defined by Zp > zi (i.e. 
the combined a and b samples) and \zp — zi\ < az, respectively, where az is the uncertainty 
in each galaxy photo-z. Samples a and b are probing the average surface density profile of 
the same set of lenses, and thus measurements of both come from the same underlying AS. 
However, because these samples have different source redshift distributions with respect to the 
lens positions, they will have different levels of lA contamination: sample b consists of lens- 
source pairs with larger Zp — zi and should thus have a lower fraction of physically associated 
galaxies. We further assume that the physically associated galaxies in both samples have 
the same average value of 71A and that the two have different levels of lA contamination 
only due to the different numbers of physically associated galaxies they contain. We discuss 
this assumption in appendix A. However, as shown in section 5.2, potential violations of it 
produce a sub-dominant bias given current levels of measurement uncertainty, and splitting 
the source sample by color should largely remove this bias. There are now only two unknown 
quantities in the lensing measurement: the true AE and the level of lA contamination per 
associated object, 71A. With two sets of lens-source pairs, we are able to solve for both 71A 
and the true AS with the lA contamination removed: 

- ( , cj-^AS^-ci^^AS, 



ci'^^cl') (AS„(i?fe - l)(Sc)i'i - AS6(i?, - l)(Sc) 



^^(^p) = m n (4-6) 

ci'')(S6-l)(Sc)S-cl'^Hi3a-l)(Se)ix) 

where we have suppressed explicit dependence on rp. 

From eq. 4.3 and the measured 71A, we can calculate the lA contribution to AS, which 
we denote ASia, and thus a fractional contamination to the observed lensing signal for a 
given source sample s. 

. . _ llK{rp){t,)^^rp){Bs{rp) - 1) 
AS AS,(rp) 

The weighting applied here for lens-source pairs is the standard scheme used for galaxy- 
galaxy lensing because it corresponds to weighting by the noise in AS. These are not the 
optimal weights for measuring the lA signal, since they give less weight to close pairs which 
are more likely to be physically associated. However, we need to solve for both 71A and AS. 
One can solve for the optimal weights to measure each of these quantities, assuming the other 
is known, and then proceed to solve for both iteratively. However, since the uncertainty in our 
signal is dominated by the AS measurement, the ideal weighting scheme should approach the 
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standard lensing weights. To test a weighting scheme that doesn't favor pairs at large hne-of- 
sight separation, we repeated our measurement using pure shape-noise weighting (removing 
the factor of and found a significantly higher uncertainty in 71A. However, as the 

precision of future measurements improves, it may be beneficial to use weights optimized for 
lA measurement with an iterative scheme. 

Note that the general technique applied here - comparing the lensing signal observed in 
two of more samples with different lens-source separations - is similar to a shear-ratio test 
to probe the geometry of the universe [59]. For the source and lens redshift distributions 
considered here, we find that varying the cosmological parameters across a reasonable range 
contributes a highly sub-dominant signal: our choice of cosmology has a negligible impact 
on I A results. However, for future studies with improved precision, the effects of assuming 
an incorrect cosmology should be considered. Conversely, if not properly treated, lA has 
the potential to contaminate shear-ratio or similar measurements and bias the resulting 
constraints on cosmological parameters. 

4.2 Estimating and reducing uncertainties 

To estimate the uncertainties in this measurement, we construct 1000 bootstrap realizations of 
all measured quantities, by random sampling with replacement, from 100 contiguous regions 
of the survey footprint. The number of bootstrap regions and realizations is chosen to ensure 
that the errors are a reasonably smooth function of scale. The quantities 71A and ASja/AS 
are then calculated for each realization, allowing us to directly determine confidence intervals 
without assuming a Gaussian distribution. Unless otherwise stated, 68% confidence intervals 
are shown. We also use these bootstrap realizations to determine model-dependent confidence 
intervals for lA. See section 5 for further discussion. 

On large scales, where S — ?• 1, the fractional error in B — 1 can become quite large. 
For source subsamples closer to the lens, where the boosts are larger, this fractional error is 
smaller. If the excess galaxies in different samples have the same composition, these boosts 
come from the same clustering. The ratio of B — 1 between these samples should then be 
constant as a function of scale, with the ratio reflecting the photo- 2; scatter into each sample. 
As discussed in appendix A, such homogeneity between subsamples is present for red galaxies, 
allowing an improvement in the boost uncertainties. We find the ratio between B — 1 for 
different samples on smaller scales (rp = 0.2-4.8/i~^Mpc), where the fractional errors are 
small. We then take the boosts for the background samples (a and b) to be the rescaled 
boosts from the assoc sample, thus providing smaller fractional errors. However, as shown 
in fig. 3, this provides only a modest improvement to the overall uncertainty in 71A. 

Furthermore, there is a systematic uncertainty in the determination of B of roughly 
3% due to large-scale structure fiuctuations and low-level variations in the lens density with 
observing conditions, which are difficult to model and are not faithfully reproduced in the 
random catalogs (see Mandelbaum et al. 2012 in prep, for details related to this lens catalog, 
and [60] for a related example in SDSS). On scales above ~ 10 /i~^Mpc, B — 1 is only a 
few percent and thus cannot be meaningfully measured. On scales between ~ 1-10 h~^Mpc, 
this systematic uncertainty dominates the errors, but the boost values are large enough for 
accurate measurement. However, this systematic uncertainty leads to some bootstrap re- 
alizations having a negative number of "excess" lens-source pairs, which is non-physical at 
these relatively small transverse separations. This artifact could impact the measured uncer- 
tainty in 7iA and fractional contamination from lA . To test the impact of this uncertainty, 
we repeat our calculations using the median boost value rather than that for each individual 
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bootstrap realization, effectively removing the systematic uncertainty, and find little effect on 
the quantities of interest. We thus conclude that the boost systematic does not significantly 
influence results on these scales. 

4.3 Identifying physically associated galaxies as "excess" galaxies 

Ideally, eq. 4.1 would be split into sums over physically associated pairs (which can have 
lA) and non-associated pairs. However, we are not able to directly measure the fraction 
of physically associated galaxies in a given sample. Instead, the boost factors determine 
the fraction of "excess" pairs. This excess is similarly described by the two-point cross- 
correlation between lens and source objects, denoted ^is- At small separations, where ^ 1, 
essentially all associated galaxies are also excess galaxies. On larger scales, however, < 1: 
the number of galaxies that would be present near the lens given a random distribution 
becomes comparable to and eventually larger than the number of excess galaxies present due 
to clustering. This complication does not reduce the applicability of the method - the physical 
mechanism of I A does not distinguish between an object that is excess above random and 
one that is not. As long as the ratio of excess galaxies to total physically associated galaxies 
is the same between the two samples, we may simply attribute the entire lA signal to the 
excess object fraction. The photometric redshift scatter should be, on average, the same for 
all physically associated galaxies, and thus this ratio should be fixed at a given Vp. While 
this approach allows for accurate lA extraction, it also can lead to seemingly counterintuitive 
results. For instance, even though alignment should be weaker at larger scales, the value of 
7iA may increase above a certain scale, as the number of excess galaxies drops more quickly 
than the total lA signal (see, e.g., [61]). The integral constraint requires that the average 
clustering across all scales be zero, while no such constraint exists for the shear correlation 
signal. This effect should be considered whenever lA is examined in the context of clustering. 

4.4 Combining correlations and photometric redshift uncertainty 

For effective comparison with both theoretical models and previous lA measurements, we 
wish to express the observables of interest in terms of underlying cosmological quantities. 
The theoretical values of both the boost factors and 71A can be written as the convolution 
of the photometric redshift scatter with the lens-source and the lens-IA cross-correlation 
functions, ^is = {6i6s) and = {6i^+). Here, 7+ = —jt, matching the standard convention 
used for shear correlation functions. 

Consider lenses at zi and sources at true redshift Zs, with photometric redshift Zp. For 
lens and source redshift distributions pi{zi) and Ps{zs) we can calculate the boost for an 
individual lens at zf. 



where H is the line-of-sight separation between the lens and source. P{zs,zi) is the weighted 
probability of an object located at Zg being placed in the source sample due to its assigned 
Zp. This quantity reflects the photometric redshift scatter into the source sample and will 
depend on the survey properties, the sources under consideration, and the definition of the 
sample: 



B{rp,zi) - 1 



J dzsPs{zs)Cis{rp,U{zs);zi)P{zs,zi) 

J dZs Ps{Zs)P{Zs,Zi) 





(4.9) 



2i 
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for a source sample with Zmin < Zp < Zmax, the boundaries of which may be defined relative 
to the lens redshift. P{zp\zs) is the probability that an object with true redshift Zg is photo- 
metrically assigned Zp. Here we have assumed that P{zp\zs) does not change as a function of 
distance from the lens due to varying average source properties or observational effects. The 
total boost factor is found by integrating over the lens distribution pi{zi): 

B(r )-l= / ^^iPijzi) I dzsPsizs)Cis{rp, U{zs); zi)P{zs, zi) 

^ f dzipi{zi) J dZsPs{Zs)P{Zs,Zl) 

In an analogous fashion, 71A can be written in terms of the and ^/g: 

- / ^ f dzipi{zi) f dzsPs{zs)^i+{rp,U;zi)P{zs,zi) wi+ 
7iA[rp) = — — TTTl , , — — r^T r ~ (^pJj 

J dzipi{zi) j dZs Ps{Zs)^ls{rp,Il; Zi)P{Zs, Zi) Wis 

where the minus sign accounts for the convention that positive "+" shear corresponds to 
negative tangential shear, and wx is the projected correlation function corresponding to (,x- 
The quantity (^J^(r))/(1 + ^is{r)) gives the lA contribution to jt per lens-source pair (both 
"excess" and "random") at r. Different lA models provide predictions for which can 
be tested against observation. As discussed in sec. 5.1, eq. 4.11 also allows for comparison 
between the results in this work and previous lA studies. 



5 Results 



We now present the results of the lA measurement technique developed above. We define 
the a and b samples using Az = 0.17 to yield a roughly equal number of lens-source pairs 
in the two and minimize measurement uiicertainty. See appendix B for a discussion of this 
split. Figure 2 shows measurements of AS for the a, b, and entire src samples.^ These 
measurements include the photo-z calibration correction factor Cz for each sample and thus 
would provide an unbiased estimate of the true AS in the absence of lA. Any statistically 
significant deviation between measurements for the a and b samples would represent either an 
error in calibration or the presence of lA. It is apparent that the measurements are consistent 
within the confidence intervals (68%), so we do not expect a statistically significant detection 
of lA. Instead, we seek to place tight constraints on the potential lA contribution. We probe 
projected separations of Vp = 0.9-10.1 /i~^Mpc. On larger scales, the systematic uncertainties 
in measured boosts, discussed above, greatly degrade the power of our method. On smaller 
scales, lensing magnification, non-weak shear, and sky systematics in the SDSS software 
pipeline limit measurement accuracy [62]. 

5.1 Previous lA measurement methods 

The work of [32] attempts to constrain lA by applying similar source photo-z cuts and then 
assuming that the background sample is effectively free of lA contamination. With this 
assumption, it is straightforward to estimate the lensing signal in the sample more closely 
associated with the lens and solve for the lA in that sample. For a background sample, b, 
with no lA: 

AS(rp) = cf)ASb(rp). (5.1) 

^Mandelbaum et al. (2012, in prep.) will include a more thorough discussion of these results, and cosmo- 
logical interpretation. 
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Figure 2. The measured AS, including the photo-z bias correction factor Cz, is shown for different 
different source photo-z samples and colors (as labeled). The agreement of the measurements for the 
different samples indicates that the effects of lA are subdominant at the current level of statistical 
uncertainty. 

Thus lA in an associated sample, a, can be expressed as 

A£a(rp) - ($) A£t{r,) 

Results using this technique are shown below in figure 3. Note, however, that direct com- 
parison with [32] is difficult. In that work, shear measurements are done directly (rather than 
using the quantity AS), and the lensing weights and photo-z bias factor are not taken into 
account. Furthermore, that work uses flux-limited SDSS Main spectroscopic sample galaxies 
[40] as lenses, rather than LRGs, which will have different lA and clustering properties. 

Other studies (e.g. [10, 11, 13]) directly measure lA using Wg+ of low-redshift SDSS 
Main spectroscopic sample galaxies as well as higher redshift LRGs (which display a much 
stronger signal). The use of spectroscopic redshifts greatly reduces the contamination from 
lensing and dilution from including non-associated pairs. Some relevant results are shown for 
comparison. It is necessary to convert such measurements to the quantity relevant for galaxy- 
galaxy lensing using eq. 4.11. This conversion is non-trivial, and thus any comparison is only 
approximate. Similarly, [63] directly measures 71A for a low-redshift sample of satellites 
around massive lens galaxies. They find a significant signal below rp = 0.1/i~^Mpc, smaller 
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than the scales considered here. On larger scales, then' results are consistent with zero lA. 
The authors of [14] constrain lA for higher-z blue galaxies {z ~ 0.6) using shapes from SDSS 
and spectra from the WiggleZ survey [64]. Their results are consistent with the limits we 
find here. 

5.2 Model-independent measurement 

Figure 3 shows the results of our lA measurement technique across a range of projected 
separations. For comparison, results are shown with and without the "extended boosts" 
technique discussed in section 4.2. Using extended boosts lowers the statistical uncertainty 
at large scales, although the difference is not significant. Also shown are the results assuming 
no lA in the background source sample (b). This measurement is more analogous to that 
made in [32], although as noted above, direct comparison with these earlier measurements is 
challenging. As expected, neglecting the lA contribution in the background sample biases the 
resulting lA measurement to a lower magnitude: the disparity between AS measurements in 
the two samples is effectively increased when allowing for lA in the background sample. The 
induced bias is smaller than the level of uncertainty in the measurement for blue, red, and 
all sources. Because this assumption also reduces the statistical uncertainty in the measure- 
ment, the resulting confidence intervals are fully contained within those of our more general 
and self-consistent measurement. Although ignoring lA in the background sample does not 
bias the lA measurement in a statistically significant way, it results in overly optimistic con- 
fidence regions. In effect, neglecting lA in the background sample has converted a statistical 
error into a systematic error. For the remainder of this work, we consider results using our 
fully consistent method. For red sources, where the assumption of a uniform excess source 
population between the two samples is justified, we use the "extended boosts," while for all 
and blue sources, we use the boosts as measured. 

Comparing measurements with and without accounting for background lA tests the im- 
pact of the assumption that 71A is the same in the two samples. As discussed in appendix A, 
the a and b samples can have somewhat different average galaxy alignment properties. How- 
ever, even if 71A differs between them, it cannot do so more strongly than in the case of 
completely neglecting lA in the background sample.^ Since the bias induced in this limiting 
case is minimal, we can be confident in applying our measurement method to all, red, and 
blue sources. 

In all cases, I A contamination is consistent with zero at the scales considered. The 
uncertainty in lA for red galaxies is smaller than for blue galaxies since the former cluster more 
strongly and have better photo-z precision, leading to a larger difference in the number of 
excess galaxies between the a and b samples. Appendix B discusses the sources of uncertainty 
and the effects of splitting sources by color and photo- 2:. 

5.3 Including minimal model dependence 

The results shown in the previous section place constraints on lA as a function of projected 
separation rp, utilizing only the lensing measurements themselves. Treating the lA mea- 
surements at different separations as independent is the maximally conservative approach. 
Including information about the underlying model of lA, which relates the signal strength on 
different transverse scales, will yield tighter constraints. Determining an accurate model for 

^In principle, the sign of the lA effect could be different in the two source subsamples, although such a 
scenario is contrived. 
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Figure 3. The measured values of 71A using three different techniques are shown (with each panel 
showing the results for a different source color sample). Black squares indicate that the original boost 
measurements have been used in eq. 4.5. Blue circles indicate the use of "extended boosts." Green 
triangles indicate that the lA contamination in the background sample was assumed to be zero. All 
three methods yield results consistent within the statistical uncertainty, although assuming that the 
background sample has no lA slightly biases the magnitude of 71A to lower values. 

lA is an ongoing theoretical challenge, made more complicated by the heterogenous nature 
of source samples used in lensing studies. For this reason, we aim to make fairly minimal 
assumptions and show the resulting constraints for two cases: a generalized power-law model 
and a model motivated from lA measurements of LRGs. 

The 1000 bootstrap realizations are combined to construct a full covariance matrix for 
the I A measurement: 

= Nk t-l ^ [7lA{ri)k7lA{rj)k - (7lA(ri))(7lA(r-j))] , (5.3) 

where angled brackets indicate an average over the realizations. We then use this covariance 
matrix to find best-fit parameters for a particular model for each bootstrap realization. At 
every value of Tp, a confidence interval is constructed using the A^'boot predictions from the 
model fits. Thus, in the case of a model with multiple parameters or one that does not 
monotonically change with its parameter, the resulting confidence region envelopes may have 
a shape different from the model itself. Because it is calculated from bootstrap realizations. 
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Figure 4. Modcl-dcpcndent confidence intervals for tiro lA signal for all, blue, and red sources are 
shown with the model- independent measurements (black data points). Our results are consistent 
with zero lA signal. Solid black (dashed green) lines denote the power-law (LRG observational) 
model. Inner lines bound the 68% confidence region while outer lines correspond to 95%. Previous 
observational results for L3 and L4 red galaxies from [10] and for WiggleZ galaxies from [14] are 
shown for comparison. These previous results have been converted to the 71A quantity, as discussed 
in section 5.2. 

the covariance matrix estimated here, C, is only an estimate of the true covariance matrix 
and will not result in a ^^-distribution (see [32]). However, since we directly construct 
confidence intervals, we only minimize and need not worry about its distribution. After 
calculating C, we make the further assumption that off-diagonal terms are zero. As discussed 
in Mandelbaum et al. (2012, in prep), shape noise dominates the covariance of the lensing 
signal on the scales used here, and thus correlations between bins become appreciable only 
when the same sources are correlated with multiple lenses - i.e., rp that is roughly twice the 
typical separation between lenses, or 20/i~^Mpc. Figure 4 shows the resulting confidence 
intervals. 

We first consider a generalized power-law model of the form 71A = Arp, where both A 
and /? are free parameters fit to each realization. This model is well-motivated in the regime 
where both Wg+ and uigg can be approximated by power laws (see, e.g., [10]). As discussed 
above, the resulting envelope of the confidence region does not follow a power-law. Instead, 
constraints are tight on the intermediate scales probed by these lensing measurements (where 
lA signal is consistent with zero), and rapidly increase on large and small scales where there 
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is little lA information. 

We also consider a model motivated by the lA observations of LRGs. Following eq. 4.11, 
7iA can be estimated as the ratio wi^/wis- Using the LRG observations of [10], we calculate 
this quantity, smoothing with a Savitsky-Golay filter and interpolating to obtain a continuous 
scale-dependence. We fit the result to our lA measurements with a single parameter for the 
amplitude, which allows for differences in lA strength and galaxy biasing between the LRGs 
and the fainter source galaxies used in this work. We use LRG observations because they 
exhibit a significantly stronger lA signal, providing sufficient S/N (above Vp = 0.9/i~^Mpc) to 
determine a well-defined scale dependence. Non-linear effects and environmental dependence 
of lA mean that LRG measurements on smaller scales are unlikely to accurately reflect the 
behavior of the relevant sources. We thus use a spline to extend to smaller scales, and note 
that constraints in this region are contingent on the smooth continuation of I A behavior. We 
fit these constraints using scales above rp = 0.44/i~^Mpc, unlike in the case of the power-law 
shape, where all measurements are used. 

This observational model assumes that the dimmer source galaxies considered here have 
the same lA scale-dependence as the LRGs, allowing for changes in amplitude due to differ- 
ences in clustering bias, alignment bias, and object ellipticity. The model is thus of limited use 
on scales where lA and clustering are qualitatively different (e.g. a different physical mech- 
anism for alignment is dominant). Unlike the power-law model, this observational model 
provides well-behaved confidence regions on larger scales because LRG observations extend 
beyond the region where we obtain constraints. As noted above, however, constraints from 
the observational model on scales below rp = 0.9/i~^Mpc are less robust. 

On larger scales, where clustering is weak, the statistical distinction between "excess" 
and "associated" pairs becomes significant, and the number of associated pairs will no longer 
be well-approximated by Wgg. In this regime, the validity of both models discussed here is 
diminished. 

5.4 Contamination to lensing signal 

As seen in eq. 4.7, the level of lA contamination in the lensing signal depends on both the 
lA signal per excess object, 71A, and the fraction of excess galaxies in the particular sample 
(directly measured by the boosts). Figure 5 shows the fractional contamination resulting from 
the lA constraints, both with and without assumptions on lA scale-dependence. Results are 
shown both for all background sources and for a subset of sources with photo-z at least 
Az = 0.17 behind the lens. Using this background sample with a photo-z margin removes a 
large fraction of excess galaxies and thus greatly reduces the potential contamination from 
lA. Due to the lensing weights used, this cut can be applied without significantly reducing 
the statistical power of the measurement. With the current level of binning in projected 
separation, the statistical uncertainty in the AE measurement is ~ 5-15%, depending on 
the sample and projected separation. For instance, at rp = l/i~^Mpc the uncertainty in AS 
is ~ 7% for all and blue sources and ~ 10% for red sources. In most cases, the limits on 
fractional contamination from lA for the scales we measure 0.1-10/i~^Mpc) are below this 
uncertainty, significantly so when model dependence is included. 

The limits on fractional contamination in figure 5 for the all and blue samples are 
conservative. As discussed in section 5.2 and appendix B, the uncertainty in 71A is larger for 
blue than for red galaxies, even though both observational results and theoretical predictions 
(e.g. [10, 23, 61]) predict that lA for blue objects is minimal. This consideration suggests that 
more realistic constraints can be obtained by taking 71A for blue galaxies to be less than that 
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Figure 5. Fractional contamination to the lensing signal is shown for all, blue, and red galaxies. Data 
points and lines follow the same convention as figure 4. The left column shows the contamination 
when the entire src sample {Zp > Zi) is used, while the right column shows contamination when 
only the background sources at greater separation {zp > zi + Az) arc used. As expected, including 
scale-dependent assumptions and applying the photo-z cut both provide tighter constraints. 

for red galaxies. Figure 6 shows the resulting fractional contamination for two limiting cases: 
7iA = for blue galaxies and 71A is the same for blue and red galaxies. With the photo-z 
cut, the more optimistic assumption (71A = for blue galaxies) yields a constraint of 1% 
on lA contamination for all source galaxies on the scales we measure. The more pessimistic 
assumption of 71A being the same for red and blue galaxies results in lA constraints of a few 
percent on these scales. 

6 Discussion 

In this work we have developed a method for measuring the intrinsic alignment of source 
galaxies in a galaxy-galaxy lensing measurement, as well as the resulting contamination 
to the lensing signal. Unlike previous lA measurement techniques, the method is fully self- 
consistent, yielding an unbiased measurement of both lA and AS, and can be used to directly 
infer contamination to the galaxy-galaxy lensing signal. Applying this method to SDSS 
LRG lenses and photometric sources, we find a signal consistent with zero lA for all, blue, 
and red source galaxies at transverse separations of ~ 0.1-10/i~^Mpc. To obtain tighter 
constraints on lA, we have assumed a functional form for its scale-dependence, using two 
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Figure 6. Fractional contamination to tlic lensing signal is shown making two different assumptions 
on 7iA. The top row shows the fractional contamination to all sources if 71A = for blue galaxies. 
The middle and bottom rows show the fractional contamination to all and blue sources, respectively, if 
7iA for all galaxies is the same as for red galaxies. These reasonable assumptions significantly tighten 
the constraints on contamination from lA. Columns and line conventions are the same as in figure 5. 

different physically reasonable models: a general power law and a model motivated by LRG 
observations. These two models yield similar results on intermediate scales, diverging on 
large scales where there are no measurements to provide information. We note that while 
our results are only directly applicable to galaxy-galaxy lensing, the alignment mechanisms 
being probed are the same as those that can contaminate shear-shear correlations. Given 
an lA model or parametrization, one could use galaxy-galaxy lensing results to determine 
constraints on the corresponding GI contamination in shear-shear correlations. The GI term 
is believed to be the primary lA systematic, since the intrinsic-intrinsic (II) term can be 
removed with redshift information. Our current work examines a redshift range somewhat 
below that most relevant for upcoming lensing surveys. Similar analysis with higher redshift 
lenses will allow us to place constraints on the GI term in future measurements of both 
galaxy-galaxy lensing and shear-shear correlations. 

lA for blue galaxies is less well constrained due to their weaker clustering and poorer 
photo- 2; precision. However, previous studies indicate that red galaxies should have a stronger 
lA signal. Thus, although a conservative approach allows for larger lA in blue galaxies, it is 
reasonable to take the results for red galaxies as an overall upper limit. 
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Compared to previous lA measurements, our analysis most closely probes the redshifts 
and galaxy luminosities relevant to current and future lensing measurements. Our results are 
broadly consistent with previous studies (e.g. [10, 13, 14, 32, 63]), which find insignificant 
lA signal for similar galaxies at lower redshifts (see figure 4). One possible exception is 
the red L3 sample of [10, 13], where there is a weak detection of lA at ~ 10 /i~^Mpc. 
However, there are notable differences between our current study and previous ones. We 
measure lA for typical sources at higher redshift than [10, 13, 32], and probe alignment 
in cross-correlation with LRGs rather than in auto-correlation. We also note that recent 
measurements [10-13] detect a significant autocorrelation lA signal for an LRG sample similar 
to our lens sample. This signal is not inconsistent with our findings, which constrain the 
I A for significantly less luminous galaxies around LRGs. However, this disparity indicates 
that different types of galaxies display different alignment in the same underlying density 
field (e.g. LRGs are significantly more aligned). Such a result is not surprising given the 
likelihood of environmental and nonlinear effects: alignment on these scales depends on more 
than simply the density field. Comparison with previous lA measurements is difficult since 
it requires assumptions on how lA and clustering behaviors change on nonlinear scales with 
object type, redshift, and environment. In addition to effects discussed above, if the primary 
alignment mechanism acts during periods of formation or accretion, followed by decreasing 
alignment due to stochastic nonlinear astrophysics, we would also expect an I A signal with 
significant redshift evolution (see, e.g., [26]). Furthermore, measuring lA at low redshift 
can be done using spectroscopic galaxies directly without having to account for the lensing 
signal, as is done in [10]. Any tension between the constraints found here and previous low-z 
measurements is therefore relevant to lA modeling but not worrisome from a consistency 
standpoint. We expect future measurements with higher S/N to detect a non-zero lA signal 
for galaxies typically used as lensing sources, although it will be weaker than that for LRGs. 

We also compare our constraints with those found in [14], who use spectroscopic red- 
shifts to directly constrain the lA of galaxies in the WiggleZ sample. These galaxies have 
a similar luminosity and roughly comparable redshift distribution to those in this study, al- 
though the WiggleZ galaxies peak at larger z. The biggest difference between the samples 
is color: WiggleZ selects emission-line galaxies with UV observations, resulting in a sample 
that corresponds to the 10-20% bluest galaxies of the "blue" galaxies examined here. Nev- 
ertheless, comparison between the results is informative. As seen in figure 4, despite the fact 
that our method is less efficient with blue galaxies than with red galaxies, the constraints we 
find for this sub-sample are tighter than the limits from [14]. Even allowing for uncertainty in 
converting between the different types of measurements and different object properties, our 
photometric approach to lA is competitive with spectroscopic methods, which necessarily use 
smaller galaxy samples. Given the wealth of upcoming photometric surveys, the disparity in 
size between photometric and spectroscopic samples with reliable shape measurements will 
only increase. 

We have placed constraints on the potential contamination to the lensing signal from 
lA: including minimal model assumptions on scale-dependence, < 5% contamination is found 
at the 95% confidence level on the scales we probe. lA thus remains a subdominant source 
of uncertainty when compared with the current statistical error on the galaxy-galaxy lensing 
signal of w 5-15%, given the binning in transverse separation. These results apply to the 
source selection and redshift quality of this particular study - contamination will be worse if 
less redshift information is used. Constraints on contamination depend not only on source 
properties such as color, but also on the photo-z selection cut used. As seen in figure 5, 
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constraints are significantly tighter when sources are selected with a photo-2: cut to remove 
a large fraction of galaxies physically associated with the lens. In this work, we apply the 
cut Zp > zi + Az, for Az = 0.17, although optimizing this selection will depend on the 
specifics of a particular survey and science goals. Applying such a cut can greatly reduce 
the possible lA contamination without sacrificing much signal, since distant lens-source pairs 
dominate the signal due to the weights. Statistical uncertainties in the AS measurement 
with and without the photo-z cut agree to within 10%. Similarly, the relative importance 
of I A is greater for lens galaxies at higher redshifts: the distribution of background source 
galaxies will be concentrated closer to the lens positions, increasing the number of physically 
associated galaxies that receive appreciable weight. 

Our findings do not suggest a clear method to cut sources by color to reduce potential 
lA contamination. Although the constraints found here are tighter for red galaxies, this is 
because our method performs better for sources that have stronger clustering and better 
photo- 2; precision. It is likely that blue galaxies display weaker alignment, and it is unclear 
to what extent this trend is offset by poorer photo-2 precision. Future lensing measurements 
with lower statistical uncertainty should resolve this issue. However, it is likely beneficial to 
remove sources with particularly high photo-z uncertainty in order to decrease both potential 
lA contamination and systematic uncertainty from imperfect photo- 2; bias correction. In this 
work, we split the source galaxies into only two sub-samples. With improved statistics, it will 
be possible to use additional sub-samples in order to break a possible degeneracy between lA 
strength and correlation between the lA properties and photometric redshift uncertainties of 
source galaxies. This issue is discussed in appendix A. 

High-precision photometric lensing science will be a primary tool in the future of ob- 
servational cosmology. As such, upcoming surveys such as KIDS, DES, HSC, and LSST 
are all designed to obtain photometry and high-resolution imaging for a large number of 
sources. While the data sets are being built and understood, early focus will likely be on 
galaxy-galaxy lensing, utilizing the large number of massive foreground galaxies with well- 
determined redshifts, such as galaxies from current spectroscopic surveys or galaxy clusters. 
In the near term, galaxy-galaxy lensing analysis is currently ongoing with the BOSS CMASS 
sample [65, 66] as lenses, and we hope to apply our method as part of the analysis. This 
sample is at higher redshift (0.4 < 2; < 0.7), allowing a probe of lA in a new redshift regime 
relevant for the next generation of lensing surveys and where lens-source overlap may increase 
the impact of lA. Similarly, we intend to use our method with cluster lenses to study lA of 
galaxies around cluster centers. Since clusters are even more highly biased than LRGs, we 
would expect larger numbers of physically associated galaxies in the source sample as well 
as a higher amplitude for alignment effects. I A contamination for cluster lenses could bias 
the measured lensing masses and the resulting cosmological results, particularly given that 
cluster lensing masses are typically determined using lensing signals on even smaller scales 
than were considered here, where lA effects will be more important. 

In the longer term, future studies should obtain galaxy-galaxy lensing signals with 
uncertainty at the ~ 1% level. Thus, although lA may become a significant systematic as 
the statistical errors decrease, applying the method outlined here will yield correspondingly 
stronger limits on lA (roughly equal to the precision of the lensing measurement) and will 
allow the removal of this contamination. When lA constraints are improved by assuming a 
specific radial scaling of the lA signal, they should reach the sub-percent level. Indeed, given 
the detection of lA in previous spectroscopic studies, these upcoming photometric surveys 
should achieve a significant detection of I A with the joint lensing and I A measurement method 
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introduced in this paper. These measurements wih further reduce systematic uncertainties 
in lensing studies. Moreover, they will probe the physical mechanisms of intrinsic alignment 
and thus improve our understanding of galaxy formation and evolution. 
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A Equating I A in different samples 

In principle, solving simultaneously for AS and 71A requires that the average value of 71A in 
the two source sub-samples, split by photo- 2; in relation to the lens position, be the same. As 
discussed above, we find that potential violations of this assumption lead to sub-dominant 
bias effects given the current level of measurement uncertainty. However, upcoming studies 
will provide greatly improved precision, and thus a discussion of this issue is warranted. 

The above assumption need not hold if some set of galaxy properties, such as color or 
luminosity, correlates with both the photo- 2; uncertainty and the level of lA. Source samples 
defined by different lens-source line-of-sight separation could then have different average 
lA properties. Since an important step in the determination of an object's photometric 
redshift is fitting to an expected spectral template based on multi-band photometry, different 
morphological types can have with different photo-z precision (see [29]). Both observational 
and theoretical studies have suggested a divide in lA properties along morphological lines: 
late-type spirals (typically blue and star-forming) and early-type ellipticals (typically red) 
are likely subject to differing physical processes that affect alignment, e.g. [23]. In the 
case of spirals, angular momentum provides the major source of object support and can 
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thus significantly influence orientation via tidal torquing. Elliptical galaxies are pressure- 
supported through velocity dispersion and are thus expected to align more closely with the 
surrounding halo and underlying tidal field. Indeed, several observations of lA have found 
much stronger signal for red galaxies ([10, 61]). 

To mitigate this issue, one can split the sample by photometric template type. Ideally, 
multiple splits would be used to isolate individual template types. However, dividing the 
source sample reduces statistical power, which quickly becomes the limiting factor. Previous 
I A studies have relied on a simple division between red and blue galaxies, and future lensing 
measurements are unlikely to consider more complicated template separation. Thus to both 
maximize the signal and make the results directly relevant, in this work we have applied a 
single split. The ZEBRA pipeline interpolates between six primary template types, allowing 
five subdivisions between each. We define all galaxies assigned to the first five interpolated 
templates, between the elliptical (Ell) and Sbc primary templates, to be red. All other 
galaxies are blue (see [29] for further information on the templates). 

The sums over weights that determine the boosts can also be used to calculate the 
effective fraction of red (or blue) galaxies in each source sample as a function of Vp, automat- 
ically accounting for the lensing weights. Figure 7 shows the effective red fraction of galaxies 
around random lenses (fr), total galaxies around real lenses (ft), and excess galaxies around 
real lenses (/e). This fraction is defined as X^red ^i/ Z^aii where the sums are taken over 
the set of random, total, or excess lens-source pairs (determined statistically). Because it 
reflects the underlying distribution and photo-z scatter of sources, with no influence from 
clustering, should be constant as a function of and is included as a reference and check 
on systematics. The quantity fe shows why splitting by color may be necessary. Since 71A in 
the two different source samples should be the same, the excess galaxies in the two samples 
should have the same properties. The excess object populations in samples a and b are quite 
different and thus will have different I A signal unless all galaxies, regardless of color, have 
the same lA behavior. 

The behavior of ft and fe seen in figure 7 is easily understood in terms of basic source 
clustering properties and photo-z uncertainties. As with many photo-z estimation proce- 
dures, ZEBRA provides more accurate redshifts for red galaxies. Poorer photo-2 precision 
for blue objects results in a larger fraction of excess galaxies (actually located at zi) being 
scattered into the b sample. This effect increases the values of fr and ft for the a sample, 
since blue galaxies have been preferentially removed. Furthermore, the fact that red galaxies 
are more highly biased than blue galaxies yields the observed scale dependence: fr and ft 
increase on small scales where clustering is more significant. The effects of this photo-z trend 
can also be seen in figure 1. The separation between Ba and Bh is larger for red galaxies, 
indicating a more accurate physical distinction between the two sub-samples. Similarly, Bi, 
for blue galaxies is larger than for red galaxies, despite having weaker clustering, because 
there are more physically associated galaxies with large photo-z error. 

To examine the effectiveness of the split, we perform an additional test using the ratio 
of boost factors between the two samples: {Baivp) — l)/{Bb{rp) — 1). The quantity {B — 1) 
measures the fractional contamination from excess sources and is given by a projection of the 
lens-source cross-correlation (eq. 4.8). Comparison of the observed boost with this prediction 
could indicate that the "excess" sources in samples a and b are statistically equivalent. In 
principle, knowledge of the source photometric and spectroscopic redshift distributions (e.g. 
from a calibration sample) allows modeling of the projection weighting {P{zs, z{) in eq. 4.8), 
from which it is possible to de-project the measured boost into a cross-correlation function. 
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Figure 7. The effective fraction of red galaxies is shown for four different source subsamples. The 
lines (nearly constant) show fr for src (solid), a (dotted), and b (dashed) samples. Open symbols 
denote /t, while filled symbols denote fe- Errors in fr and ft are negligible and are thus not shown. 
The composition of the a and b samples are clearly different. 



However, it is impractical to obtain P{zs, zi) with adequate resolution to probe the correlation 
function on relevant scales {vp < 10 /i~^Mpc). Furthermore, it is challenging to construct a 
set of galaxies from the calibration set that reliably measures the photo-z behavior of sources 
near lens galaxies. Another way to probe the sufficiency of the object split is to look at the 
ratio of measured B — 1 values for different samples. Because of large photo-z uncertainties, 
P(zs,zi) will be broad compared with and is thus effectively constant in the numerator 
of eq. 4.10. If the and dependencies of S,is separate, we find: 

^^"'^ - B^irp) - 1 ^ w^s^ir^y ^^-^^ 

where wis{rp) = / dnps(n)^is(rp, 11) is the projected cross-correlation function, and Sa and 
Sb denote the sources found in samples a and b. The constant of proportionality involves 
the ratio of integrals over lens and source distributions as well as the photo-z projection 
weighting P. 

If samples a and b contain a statistically similar collection of sources, this ratio will 
be roughly constant across the range of scales we wish to probe. Scale-dependence in the 
ratio could indicate that the samples have a different composition and thus have different 
scale-dependent biases on the cross correlation. 

In figure 8, the quantity R{rp) is plotted, for all, red, and blue source galaxies. Previous 
work indicates that red galaxies have reasonably well-understood patterns of small-scale clus- 
tering, photo-z error, and lA [13, 67, 68]. The roughly constant value of R{rp) measured for 
red galaxies supports the fact that we have selected a source sample with uniform properties 
across the photo-z split. The lA properties of "blue" galaxies are poorly understand, allowing 
for more significant variation within the broad sample. The measurement of R{rp) for blue 
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Figure 8. R{rp) is shown for different object colors (as labeled). The ratio is constant, within the 
estimated uncertainty, for all and red sources but not for blue sources. 



sources in figure 8 indicates that the sample exhibits significant heterogeneity in clustering 
(and hkely lA) properties. Results for blue sources are thus more subject to systematic un- 
certainties, although as discussed above, these uncertainties are currently sub-dominant to 
the relatively large statistical uncertainty on intrinsic alignment contamination. 

However the R{rp) diagnostic is not conclusive. For instance, the entire source sample 
with no color split contains multiple galaxy populations, each with different clustering and 
lA properties. The value of R{rp) depends on the relative abundance as well as the scale- 
dependent bias of each. As seen in figure 8, these factors result in a roughly constant value 
of R{rp) for the entire source sample, despite its composite nature. Similarly, if assumptions 
on the relevant scale dependencies in the measured boost are violated, a uniform sample 
could exhibit a non-constant Rivp). We reiterate that for this work, even if the assumption 
of sample homogeneity is violated, as we expect it to be with the all and blue samples, 
our method still provides a more self-consistent and unbiased lA measurement than earlier 
techniques. The bias potentially introduced by the assumption is shown in section 5.2 to 
be small. In future studies with higher measurement precision, the implications of this 
assumption should be considered. 

More generally, with sufficient statistical power, it will be possible to determine the effect 
of correlations between lA properties and photometric redshift uncertainties. We expect such 
correlations to exist, even within a single morphological category, since both lA and photo-z 
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quality are affected by galaxy luminosity (e.g. [10, 55]). Splitting the source galaxies into 
more than two sub-samples will allow a measurement of how 71A varies. Any two sub-samples 
can be used to solve for 71A, and the results from different sub-samples can then be compared. 



B Sources of uncertainty 



The sources of uncertainty in the lA measurement are seen in eq. 4.5. AS is subject to shape 
and measurement noise. Scatter in the number of random-source pairs leads to uncertainty 
in both the boost factors and (Sc)cx- All of these quantities are affected by photo- 2; scatter, 
the bias from which is corrected with the Cz factor. The uncertainty in Cz, as measured 
from calibration sets, is at the ~ 2-3% level [29], including both statistical and systematic 
effects, and we thus do not include it. We note, however, that 71A depends only on the ratio 
of Cz between the source samples. Thus, for deeper surveys for which it can be difficult to 
obtain sufficiently representative photo-z calibration samples, the uncertainty in Cz can be 
mitigated by measuring the ratio of the lensing signal for the two samples on large scales 
(~ 50/i~^Mpc), where the only difference should be due to photo-2 calibration. With the 
exception of c^, all the sources of uncertainty mentioned here are automatically included in 
the bootstrap realizations. On the range of scales considered here, the most significant source 
of error is the shape and measurement noise in AE. 

As seen in figure 2, AS is measured at similar levels of precision for red and blue 
sources, since they have roughly the same number of lens-source pairs and measurement 
noise per pair. However, the uncertainty in the lA signal for blue galaxies is higher than for 
red galaxies. This difference is due to stronger clustering and more precise ZEBRA photo-z 
measurement for red source galaxies. The denominator of eq. 4.5 converts the difference 
between AS measurements to the lA signal. When Ba ^ B^, the denominator scales as 
Ba — 1, and uncertainty is reduced when strong clustering yields a large value of Ba- If the 
difference between Ba and Bi, is smaller, as would be the case for large photo-z uncertainties, 
the denominator is further reduced. These two factors combine to degrade the precision of 
the lA measurement for blue galaxies. 

Examining these sources of uncertainty also allows us to understand the effects of differ- 
ent source splitting on the measurement precision. If boost factors were held fixed and were 
perfectly known, the error in ASq — AS^ would be minimized when the a and b samples 
had equal lens-source pair numbers, weighted by (cz-B(Sc))^^ for each sample. However, 
changing the location of the photo-z split will affect both the number of lens-source pairs 
and the overall boost in each sub-sample. Increasing the value of Ba will decrease the over- 
all error for a given measurement uncertainty. Moreover, uncertainty in the boost factors 
themselves contribute a non-negligible error, which decreases with larger Ba in the relevant 
regime. Finally, having (Ba — 1) ^ {Bi, — 1) both decreases the uncertainty and avoids 
a singular lA estimator, which would yield a large, non-Gaussian variance. Our choice of 
Az = 0.17 was motivated by these considerations, although it is not rigorously shown to 
be the optimal value. We tried other, similar splitting schemes, and found our choice yields 
smaller uncertainties in 7ia. 
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